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We consider two empirical relations using data only from the prompt emission of 
^ . Gamma- Ray Bursts (GRBs), the peak energy {Ep) - peak luminosity {Lp) relation (so 

^ ' called Yonetoku relation) and the i^p-isotropic energy (iSiso) relation (so called Am- 

\ ati relation). Both relations show high correlation degree, but they also have larger 

dispersion around the best fit function rather than the statistical uncertainty. Then 
we first investigated the correlation between the residuals of Lp and Eiso from the 
best function, and found that a partial linear correlation degree is quite small of 
' Pip Eisa-Ej, = 0.379. This fact indicates that some kinds of independence may exist 

I between Amati and Yonetoku relation even if they are characterized by the same 

QQ , physical quantity Ep, and similar quantities Lp and Ei^o which mean the bright- 

ness of the prompt emission. Therefore we may have to recognize two relations as 
the independent distance indicators. From this point of view, we compare constraints 
on cosmological parameters, r2„i and Q\, using the Yonetoku and the Amati rela- 
' tion calibrated by low-redshift GRBs with z < 1.8. We found that they are different 

, in 1-CT level, although they are still consistent in 2-a level. In this paper, we intro- 

duce a luminosity time defined by = Eiso/Lp as a hidden parameter to cor- 
rect the large dispersion of the Yonetoku relation. A new relation is described as 
(Lp/1052 erg s-i) = 10-3-87±o.i9(£;p/keV)i-82±o o8(Ti/s)-0-34±o.09 We succeeded in 

' reducing the systematic error about 40% level, and might be regarded as "Funda- 

I mental plane" of GRBs. We show a possible radiation model for this new relation. 

Finally, applying the new relation to high-redshift GRBs with 1.8 < z < 5.6, we 
obtain (Q^rm^h) = (O-IT^qqs' ^-^l-o 6i)' which is consistent with the concordance 
cosmological model. 
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1 INTRODUCTION 

In ou r previous papers l|Kodama et al.l l2008l : iTsutsui et al.l 
I2OO9I ). we calibrated the relation between peak energy Ep 
and peak luminosi ty Lp of prompt GRB emission (so called 
Yonetoku relation (jVonetoku et a l. 2004)) by 33 low-redshift 
GRBs {z < 1.62) whose luminosity distances were es- 
tima.ted from SNela llRiess et al.l l2007l : IWood-Vasev et al.l 
I2OO7I : iDavis et al1l2007h . Then we used the calibrated Yo- 
netoku relation as a distance indicator like the period- 
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luminosity relation of Cepheid variables, and exten ded the 
Hubble diagram up to z = 5.6. lTsutsui et al.l l|2009l ) showed 
that GRBs constrain cosmological parameters in a different 
way from SNela, and GRBs could be useful to probe cos- 
mological expansion of high-redshift universe where no SNIa 
has been observed. 

The Amati relation l|Amati et al.ll200^ is another rela- 
tion for prompt emission property. It involves peak energy 
Ep and isotropic energy -Eiso, and was originally derived un- 
der a given set of cosmological parameters. Therefore, the 
circularity problem arises if one applies naively the Amati 
relation to determine cosmological parameters. To overcome 
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this difficulty, in this paper, we first cahbrate the Amati re- 
lation as we did for Yonetoku relation without assuming any 
cosmological models but using luminosity distance given by 
SNela for z < 1.8. 

Although the strong correlation between Eiso and Lp 
are generally confirmed, we found little trend between their 
data residuals of ALp and Ai7iso from the best fit function of 
each relation. We suggest statistical independence of Amati 
and Yonetoku relations while these two relations have the 
same parameter Ep (§-2). 

Next, we will extend the Hubble diagram with the ob- 
tained relations and make constraints on density parame- 
ters, (f2m,f^A). It is shown that the two Hubble diagrams 
differ systematically at high redshifts and, as a result, two 
different constraints are obtained (§-3). Although the differ- 
ence is not so significant (l-cr level), taking the relatively 
large systematic errors in the two relations themselves into 
consideration, it is suggested that there may be a hidden 
parameter which characterizes the prompt emission and re- 
duces the systematic error of the distance indicator. We in- 
troduce the luminosity time {Tl = Eiso/Lp) as a possible 
hidden parameter to derive the Ep-Th-Lp relation, and put 
constraints on cosmological parameters (§-4). Finally we give 
some comments on other relations used as distance indica- 
tors in the past, and argue the advantage of our new relation 
(§-5). Throughout the paper, we fix the current Hubble pa- 
rameter as Ho = 66 km s~'^Mpc~'^. 

We note that, in whole this paper, tentatively we do 
not take possible selection effects and evoluti on effects on 
relati ons into accounts although a hot debate jButler et al.l 
I2OO7I : Q I2OO7I : iBasilakos fc Perivolaropouloj 120081 ) exists^ 
We think that the final resolution needs the increase of the 
events with enough well determined parameters such as red- 
shifts, Ep, Lp and Eiso since we should divide all data into 
several sub groups to check possible selection effects and 
evolution effects. 



2 CALIBRATION OF AMATI RELATION AT 
LOW REDSHIFTS 

The typical spectrum of the prompt emission of GRBs can 
be expressed as expone ntially connected broken power-law, 
so called Band function ijBand et al.l ll993^) . Then we can de- 
termine spectral peak energy Ep, corresponding to the pho- 
ton energy at maximum in vF,j spectra. There are two em- 
pirical relations that relate prompt emission pro perty with 
Ep. E p-Eiao relation is the first one found by Amati et al.l 
(|2002l ). which connects Ep with the isotropic equivalent 
energy -Eiao. The seco nd one is Ep-Lp relation found by 
lYonetoku etaP (l2004h wh ich was used in our previous pa- 
pers (|Kodama et al.ll20ol : iTsutsui et aLll2009l ) . 

We first calibrate Amati relation in the same way as 
in ou r previous papers ijKodama et al.l [20081 : pTsutsui et al.l 
I2OO9I ) and analyze the correlation of the residuals of GRB 
data from the relations and the partial correlation coeffi- 
cient. If there are a small degree of correlation between Lp 
and Eiso after removing the effect of Ep, it suggest the possi- 
ble independence of the distance indicators, and might sug- 
gest the existence of a hidden parameter common to each 
relation. 

We found an empirical formula for the luminosity dis- 



tance as a function of redshift from 192 SNela observa- 
tions (|Riess et al.ll2007l : lDavis et al.ll2007l : IWood-Vasev et all 
|2003), 
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= 6.96 X z^-''^ + 14.79 x 



(1) 



The reduced chi-square of the formula is ~ 0.995. Note 
here that the formula is not unique and a different formula is 
possible. Note also that we do not assume any cosmological 
models at this stage, but simply assume that the Type la 
supernovae are the standard candles for 0.168 < z < 1.755. 
Furthermore, we neglect the errors in Eq. ([TJ in the follow- 
ing analysis, which leads to the underestimation of errors 
in cosmological parameters. Our purpose here is to compare 
distance indicators and find a better indicator so that this 
neglect would be reasonable. 

We apply this formula to 31 lo w redshifts GRBs with in 
z < 1.6. ( For details of our data, see lKodama et all l|2008l ).) 
In Fig. [T] we show the peak energy Ep and the isotropic 
energy _Eiao of 31 GRBs with z < 1.62. The solid line is the 
calibrated Amati relation given by. 
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where the statistical errors are indicated and Eiso = 
47rci2 S'boi/(l + 2) where Sboi is the bolometric fluence es- 
timated in 1-10000 keV energy range in GRB rest frame. 
The Pearson correlation coefficient is 0.943 and the re- 
duced chi-square is ~ 28.5/29 with the systematic 
error asys — 0.35. Here we include not only errors in 
Eiso but also errors in Ep, and the systematic error of 
this relation (Jsys, so the chi-square function is defined 
as x'(^,B) = S(logS°i'<f - A BlogEpf/{al,,s + 
where the weighting factor al^^as = (^logEi,^ + B^^^iogEp- 
The value of asys can be estimated by the value such that a 

fit to the Amati relation calibration produces a reduced 

of unity. 

This kind o f systematic error wa s also found in the Yo- 
netoku relation (|Kodama et al.ir2008l ). 
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where Lp is 1-second peak luminosity. This Yonetoku re- 
lation is slightly different from that in the previous work 
because we include not only Lp error but also Ep error, and 
the systematic error in this paper. The Pearson correlation 
coefficient is 0.948, the reduced chi-square is x^ = 30.6/31, 
and the systematic error is asys = 0.27. The Amati relation 
has slightly larger systematic error than the Yonetoku rela- 
tion and, in both cases, systematic errors are significantly 
larger than measurement errors. 

Fig.[2]shows the correlation between two residuals which 
are defined as ALp = (logL"^^ ~ logL^^P) and A£;iso = 
(log S,°^=-log B^^^P). Here, Lf^ {E?^;) is the observed value, 
and L'^p^ (Eiso) is the expected quantity from Yonetoku and 
Amati relations for observed Ep, respectively. We can see 
that there is no correlation between ALp and AEiso- We 
performed a linear partial correlation test for these 31 sam- 
ples, and found the partial correlation coefficient is quite 
small of pLp Ei^„-Ep = 0.379. Here P12 3 means the correla- 
tion coefficient between the first and the second parameters 
after fixing the third parameter which is controlling two cor- 
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Figure 1. The peak energy {Ep) and isotropic energy (Biso) of 
31 GRBs with z < 1.62. The solid line is the calibrated Amati 
relation given by Eq l[2]l and dashed lines represent the 1-sigma 
region. The systematic error of this relation is reduced to asys = 
0.35 which is larger than measurement errors. 
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Figure 2. Residuals from calibrated Yonetoku relation and Am- 
ati relation from the observed Ep. Here, the residuals are the dif- 
ferences between the observed quantities, L'^^ and E°^^ , and the 
expected quantities, and Ef^^, from Yonetoku and Amati 

relations for the given observed Ep, respectively. There seems to 
be no correlation, which implies the two relations are independent 
distance indicators. 



relations (the parameter 1-3 and 2-3 relations). This fact 
indicates that two distance indicators may be independent 
from each other. In the next section (§-3), we investigate 
whether the cosmological parameters measured with two in- 
dependent relations are consistent or not. After that, in §-4, 
we explore a hidden parameter which is the cause of the 
intrinsic dispersion of Amati and Yonetoku relations. 



i 45 - 




2 3 4 

redshift 



Figure 3. Extended Hubble diagram from Yonetoku relation 
(blue) and Amati relation (red). A systematic difference seems 
to exist in high redshift, although it doesn't seem in low redshift 
GRBs. 



3 CONSTRAINTS ON COSMOLOGICAL 
PARAMETERS 

We apply these Amati and Yonetoku relations to 29 GRBs 
with high redshifts, 1.9 < z < 5.6, to determine the luminos- 
ity distance as a function of z. We include not only measure- 
ment errors of Eiso (Lp), and Ep but also systematic errors 
we estimated in §-2 and the 1-sigma uncert ainties in Eg. ^ 
(Eq. Q) to estimate luminosity distance (jSchaeferl l2007l l . 
For example, in case of Amati relation, we define y,Q{zi) and 



A*o(2i) 



as below, 

5 
2 

25 

4 



[A + B\ogEp-\og[^ 



(Ja + {o'B log Ep)'^ + a. 



iTvSbol 



+ z 

2 



(4) 
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Fig.|3]shows an extended Hubble diagram up to 2: = 5.6 
from Amati relation (red) and Yonetoku relation (blue). A 
systematic difference between red and green points seems to 
exist especially in high-redshift region. 

Then we derive constraints on cosmological parameters. 
In the A-CDM model with flk = fim + ^^A — 1, the luminosity 
distance is given by, 

d*(z, fl,n, Oa) 

: sm{y^F{z)) if fifc > 



(1 + ^) 



: sinhi^/^IhEiz)) if < (B) 
if Qk=0 



with 

F{z) 



dz 



^■m{l+z'f -ilk{l + z'Y 



-1/2 



The likelihood contour is defined by. 



^o(zi) — ^J^'^{zi, fim, ^A, ) 



2 

Xbcst ? 



(7) 



(8) 



where /^"'(zi, fim, ^a) = 5 log(dl''/Mpc) + 25 and xLst rep- 
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Figure 4. Constraints on (Hmif^A) plane from Amati relation 
(red) and Yonetoku relation (blue). The contours correspond to 
68.3% confidence regions, and black solid line represents the flat 
universe. They are slightly different, although they are consistent 
in 2-cr level. See also Table [T] 
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n 25+0-31 
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9.36/28 
9.71/29 


Ep-TL-Lp (flat) 


17+0-15 
24+0-11 


-1 9-1+0.07 
^■^^-0.61 


16.59/27 
17.50/28 



Table 1 . Constraints on (f2m , A ) in non-fiat and fiat universe 
with 1-(T errors from Amati, Yonetoku and Ep-T^-Lp relations 
and their reduced chi squares. Constraints from Amati and Yo- 
netoku relations are inconsistent in l-cr level. 



resents the chi-square value for the best-fit parameter set of 
Q.m and 

In Fig. |4j we show the likelihood contour from Amati 
(red) and Yonetoku relations (blue), and the best-fit values 
with 1-a errors are shown in Table [T] Interestingly, they are 
slightly different, although they are consistent in 2-a level. 



eters, because the systematic errors in the Amati and Yo- 
netoku relations are larger than their mea surement errors. 
(|Ghirlanda et al.ll2006l : iFirmani et "al]|2006l ) There are some 
studi es using jet break tim e in afterglow as such a param- 
eter l|Ghirlanda et al.ll2004l : iLiang fc Zhand 120051 ). but the 
number of GRBs which have observed jet is small. Here we 
restrict our discussion in a prompt emission property. In the 
past studies, prompt emission was characterized by a time 
scale, the duration o f most intense parts of the GRB (ro.45) 
JPirmani et al.|[2006l : iRossi et"al 1 l2008l :[ CoUazzi fc Schaeferl 
l20Q8t ). Here w e adopt a time sc a le call ed luminosity time Tl 
introduced bv lWiUingale et all l|2007t) as. 



Tl = 



(9) 



{l + z)Fp 

The luminosity time does not depend on detector's energy 
band because it is defined by bolometric flux and fluence. 

We assume the correlation among Ep, and Lp to 
be of the form, log Lp = A-\- BlogEp -\- ClogT^. Then we 
obtain. 
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from low redshift 30 GRBs in 0.16 < z < 1.7. In Fig. [5] 
we show Ep-TL-Lp relation. The Pearson correlation 
coefficient is 0.971, the reduced chi-square is = 26.9/27 
with the systematic error u^^" = 0.15. We include 
not only errors in Lp but also errors in Ep, and Tl 
so the chi-square function is defined as 'x^{A,B,C) = 



E(log LI"' - A- BlogEp^C log Ti,)7(a, 



2 

meas 



+ 



s) 



where the weighting factor a„ 



(1 + 2C)a' 



{BaiagEpY + (C'^iogTi)^ , and aiogTZ, is estimated by using 
the error propagation equation without a crossterm between 
Lp and Eiso- The factor 2C in the front of cfogLp comes 
from the fact that the definition of Tl includes Lp. Note 
that the contribution of the additional error term ot^ to 
chi-square value is little because of its small slope. Thus, we 
conclude that the additional term really improves relations. 
The systematic error is substantially reduced compared to 
those of Amati and Yonetoku relations, and now compa- 
rable to the measurement error. Thus this relation could 
be re garded as "Fundamental plane" (jPio rgovski fc Davis! 
Il987f) of GRB prompt emission. 

Here we excluded one outlier, GRB070521, from 
the fitting of Eq. (|10|) . Actually, the host galaxy of 
G RB070521 is detected inside an er r or cir cle of XRT 
by iHattori. T.. Aoki. K.. fc Kawai, nTI (|2007l ) using Sub- 
aru Telescope after 40 minutes from the trigger, but they 
couldn't detect bright afterglow. Thus, the real redshift 
may be larger, which is why we exclude GRB070521 from 
our analysis. Note that the probability of the miss iden- 
tification of the host galaxy fr om only XRT observation is 
about 7% (|Cobb fc Bailvnll2008l ). Futhermore there might be 
another population of GRBs so that the secure classifica- 
tion of GRBs is one of the current issue of GRB's study. 



4 NEW RELATION 

In this section, we seek for a hidden parameter which re- 
duces the systematic error in the relation. Note that the 
discussion in this section does not concern with indepen- 
dence of the two relations. If we could find this parame- 
ter, we can make strong constraint on cosmological param- 



iTsutsui et al. I l|in pretj ). 

We show a possible derivation of the new relation 
(Eq. (|10|l under t he photosp heric model of the prompt emis- 
sion of GRBs ( lloka et al] (|2007i ) and references therein). 
The luminosity is given by 

(11) 
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in 31 GRBs with z < 1.62. 



Figure 5. The Ep/r0 i87 and Lp 
The correlation is improved than Amati and Yonetoku relations. 
The solid line shows the best-fit curve without one outlier (green 
square:GRB070521), and dashed lines represent the 1-sigma re- 
gion. The Pearson correlation coefficient is 0.971 and reduced chi- 
square is = 26.9/27 with the systematic error <Jsya = 0.15 
which is comparable to the measurement error. 
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Figure 6. Constraint on {Qm, Qa) plane from Ep-T^-Lp relation. 
See also Table [T] 



where r, F and T' are the photospheric radius, the gamma 
factor and co-m oving temp e ratur e of the photosphere, re- 
spectively. Since lloka et aP ()2007l ) assume that the energy 
is supplied by the relativistic collision of the rapid shell of 
mass rrir and Lorenz factor 7r with the slow shell of rus 
and 7s. Then under the perfectly inelastic collision model, 
r is given by = (mr7r + ms7s)/(m,./7r + m.s/'ys)- Since 
7^ ^ 7s and mr'jr + rriajs oc Eiao, we can reduce oc 
Eiaoja/ms. If we regard Ep ~ FT' , we can rewrite Ea. (|lip 
as L oc r'^Ep/r'^ oc r^nisEp/ Eisojs- Now let us assume r, 
rris and 7s 7l are constants. Then we have L oc EpT£^'^^ . 
The above relation is essentially the same as Ea. (|10p if we 
consider 2-cr error of the power index in Eq. (|10|) . 

Possible r easons for the ab ove three assumptions are 
as follows. In lloka et al.l (|2007l ) model r is similar to the 
radius of the progenitor star so that it could be constant. 
Tl can be regarded as the effective duration of the burst. 
Then cjITl is the radius that the last rapid shell catches 
the slow shell and we expect that this is also the order of 
the radius of the progenitor star, which is constant. We have 
no reason why nis is constant. However if ttIs obeys the log 
normal distribution like the other observables in GRBs we 
may regard it essentially constant. If these assumptions are 
reasonable, the new relation (Eq. (|10|l ') could be derived in 
the photospheric model of the prompt emission of GRBs like 
lloka et al.l (|2007l ). 

Finally we use this new relation to put constraints on 
cosmological parameters. Here again we include not only 
measurement errors of Lp, Ep, Tl but also systematic er- 
rors and the 1-sigm a uncertainties in Ea. (|10|l to estimate 
luminosity distance (|Schaefedl2007h . 

The concordance cosmology is still consistent in 1- 
(j level. The constraints on cosmological parameters are 
(n^,nA) = (0.17i;5:J^,1.21«0{) and xl = 16.59/27 for 
non-flat universe, ri„ = 0.241°:^ and xl = 17.50/28 for 
flat universe. (See Table [1]) 



5 DISCUSSION 

Recently some authors cxtondod Hubble di agram up to 
z ~ 6 using vari ous luminosity in dicato rs ([ Amati et al 



200? 



Liang et all [ioosi: ISchaefej l2007l: iFirmani et al 



20061: IGhirlanda et all I2OO6I). iGhirlanda et al.1 l|2006l ). 
Firmani et all l|2006l ). and ISchaefeij (|2007l ) are pioneermg 
works for GRB cosmology, but they are caught in circularity 
pro blem because th ere are few GRBs at low redshift. 

ISchaefed (120071 ) obtained cosmological constraints from 
lag (riag)- luminosity relation, variability (y)-luminosity re- 
lation, i?p-jet coUimated energy {E-^) relation (so called 
Ghirlanda relation) , minimum rise time (Tp tT)-luminosity re- 
lation, and Yonetoku relation. iLiang et al.l ifioOS ') calibrated 
these relations by luminosity distances from SNIa observa- 
tions. 

However, iTsutsui et al.l (|2008l ) showed redshift depen- 
dence of Tiag-Lp relation analyzing 565 BASTE GRB sam- 
ples with pseudo-redshifts estimated by Yonetoku relation. 
This suggests that the relation cannot be used as a dis- 
tance indicator. The dependence might come from the fact 
that TiagS were evaluated from several fixed energy bands 
depending on the detectors and they are different energy 
in GRB rest frame ,so it sufi'er from K-correction problem. 
This argument applies to Vs and trts. 

Although Ep-E-y relation and Liang & Zhang relation 
have much smaller systematic errors than both Yonetoku 
and Amati relations, there are many GRBs without a jet 
break or with multiple jet breaks (missing or multiple jet 
break problem) so that it is not certain whether the jet 
break can be used to characterize GRB emission. In con- 
trast, our new relation does not suffer from these problems, 
because it is totally defined by the prompt emission prop- 
erty and it has as a small systematic error as Ep-E^ and 
Liang & Zhang relations. Thus we could expect that the new 
relation would put GRB cosmology to the next promising 
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stage, as Phillips relation (|Phillipslll993) and Fundamental 



plane l|Diorgovski fc Davis 



1987h did for SNela and ellip- 



tical galaxies. However, we must emphasize that we need 
detailed studies of the new relation with much larger num- 
ber of GRBs and examination of systematic errors in order 
for GRB to be regarded as a reliable tool for cosmology 
like SNIa, cosmic microwave background, baryon acoustic 
oscillation and gravitational lens. Now ongoing Missions like 
Smft , Fermi and Suzaku , and the collaboration of many 
observer on ground will promise the progression of GRB 
cosmology. 
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